NEAR-PRIMITIVE ROOTS 



PIETER MOREE 

Abstract. Given an integer t > 1, a rational number g and a prime p = l(mod t) 
we say that g is a near-primitive root of index t if i'p{g) — 0, and g is of order 
(p ~ modulo p. In the case g is not minus a square we compute the density, 
under the Generalized Riemann Hypothesis (GRH), of such primes explicitly in 
the form p[g)A, with p[g) a rational number and A the Artin constant. We 
follow in this the approach of Wagstaff, who had dealt earlier with the case where 
g is not minus a square. The outcome is in complete agreement with the recent 
determination of the density using a very different, much more algebraic, approach 
due to Hendrik Lenstra, the author and Peter Stevenhagen. 



1. Introduction 

Let g G Q\{ — 1, 0, 1}. Let p be a prime. Let Vpi^g) denote the exponent of p in the 
canonical factorization of g. If i'p{g) = 0, then we define 

rg{p) = [{Z/pZ)* : {g mod p)], 

that is rg{p) is the residual index modulo p of g. Note that rg{p) = 1 iff is a 
primitive root modulo p. For any natural number t, let Ng^t denote the set of primes 
p with i^pig) = and rg{p) = t (that is Ng^t is the set of near-primitive roots of index 
t). Let S{g, t) be the natural density of this set of primes (if it exists). For arbitrary 
real x > 0, we let Ng^t{x) denote the number of primes p in Ng^t with p < x. 

In 1927 Emil Artin conjectured that for g not equal to —1 or a square, the set 
Ag 1 is infinite and that Ag i(x) ~ CgATi{x), with Cg an explicit rational number, 



^ = n (l - ^(^) - 0-3739558, 



V 



and 7r(x) the number of primes p < x. The constant A is now called Artin's constant. 
On the basis of computer experiments by the Lehmers in 1957 Artin had to admit 
that 'The machine caught up with me' and provided a modified version of Cg. See 
e.g. Stevenhagen [12] for some of the historical details. On GRH this modified 
version was shown to be correct by Hooley [2]. 

Thus 5{g, 1) is explicitly known (under GRH). Determining similarly 6{g,t) turns 
out to be rather more difficult and for ease of exposition we first consider the case 
where g > 1 is square free. In this case work of Lenstra [3] and Murata [H] suggests 
the following conjecture (with as usual fx the Mobius function and Cfe = e^'^*/''). 
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Conjecture 1. Let g > 1 be a square free integer and t > 1 an integer. The set 
Ng^t has a natural density S{g,t) which is given in Table 1. We have 

Ng^t is fimte iff 6{g,t) = iff ^ = l(mod 4), 2 ft, g\t. 

We note that if g = l(mod 4), 2 \ t and g\t, tlien Ng^t is finite. To see this note 
that in this case we have (|) = 1 for the primes p = l(mod t) by the law of quadratic 
reciprocity and thus rg{p) must be even, contradicting the assumption 2 \ t. 
Note that if a set of primes is finite, then its natural density is zero. The converse 
is often false, but for a wide class of Artin type problems (including the one under 
consideration in this note) is true (on GRH) as first pointed out by Lenstra [3]. 

Given an integer a and a prime q, we write to denote the g-part of a (that is 
Og = q'^ with q^\a and q'^^^ \ a). We put 



B{g,t)= n 



-1 



A 



n 



p 



1 
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Note that if g\t, then in the definition of B{g,t) we have the empty product and 
hence B{g,t) = 1. It follows that if further t is odd and g = l(mod 4), then 
5{g, t) = 0. The maximal value of 6{g, t) that occurs is 2E{t). Table 1 we took from 
a paper by Murata |9j . We will show that the densities in Table 1 can be compressed 
into one equation, namely ([7]). 

Table 1: The density 6{g,t) of Ng^t (on GRH) 



9 


t2 


5{g,t) 


g = l(mod 4) 


t2 = l 
2\t2 


{l-B{g,t))E{t) 
{l + B{g,t))E{t) 


g = 2 (mod 4) 


t2<A 

t2 = 4 
t2 > 4 


E{t) 

{l-B{g,t)/3)E{t) 
{l + B{g,t))E{t) 


g = 3 (mod 4) 


t2 = l 
t2 = 2 

t2 > 4 


E{t) 

{l-B{g,t)/3)E{t) 
il + B{g,t))E{t) 



Theorem 1. Conjecture^ holds true on GRH. 
The proof is postponed untill Section |2l 

2. Generalization to rational g 

A natural next problem is to study what happens if one relaxes the condition that g 
should be square free. Our starting point here will be a result due to Wagstaff [13j . 
We need some notation. We put 

fi{n){nt, h) 



S{h,t,m) = ^ 



mint 



ntip{nt) 



with (f Euler's totient function. 
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Theorem 2. [13j. (GRH). Let g G Q\{— 1, 0, 1} and t > 1 be an arbitrary integer. 
Write g = ±gQ, where go & Q is positive and not an exact power of a rational and 
h > 1 an integer. Let d{go) denote the discriminant o/Q(y^). The natural density 
of the set Ng^t, 6{g,t), exists and is given by 



n=l 



fi{n) 



nt, 9 



1/nt] 



which equals a rational number times the Artin constant A. Write go = gig2, where 
gi is a square free integer and g2 is a rational. If g > 0, set m = lcm{2h2, d{go)) ■ 
For g < 0, define m = 2gi if2\h and gi = 3 (mod 4), or/i2 = 2 and gi = 2 (mod 4); 
let m = \cm{4:h2,d{go)) otherwise. If g > 0, we have 6{g,t) = S{h,t,l) + S{h,t,m) . 
If g < we have 



6{g,t) = S{h,t,l] 



^S{h, t, 2) + ^S{h, t, 2h2) + S{h, t, m) 



(3) 



In case g > or 2 \ h, Wagstaff expressed S{g,t) as an Euler product. By the 
work of Lenstra [3] we know this is also possible in general. The next theorem 
achieves this. Partial inspiration for it came from recent joint work with Lenstra 
and Stevenhagen, see Section El 



Theorem 3. (GRH). Let g G Q\{ — 1, 0, 1} and t > 1 be an arbitrary integer. Write 
g = ±(7o ' where go & Q is positive and not an exact power of a rational and h > 1 
an integer. Let d{go) denote the discriminant ofQ{y/go). Put Fp = QiCp^ 9^^^) ■ Put 

P 



A{g,t) 



t^ 



p\t, hp\tp 



Put 



Put 



Hi 



n 



-1 



p\d{go), p\2t 



E2{m2) 



-1/3 



[Fp : Q] - 1 

if ^21^2; 

if m2 = 2t2 7^ 2; 
if m2 = 2t2 = 2; 
if 1712 t 2t2, 



(4) 



We have 

A{g,t) it,h) 



A P^ -p 

p\th 



n^^T^ 11^(^-1) iK^'-i) n ^(^-2) 



(5) 



Ptp\hp 



p\t 
hp\tp 



p\h, pt*l 



where 



2t if ^ < 0, 2\h, 2\t; 
t otherwise. 



Note that A(g, t) = iff g > 0, 2\h and 2 \ t. 

The natural density of the set Ng^t exists, denote it by 6{g,t). 
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Put Vq = lcm(2/i2, d{gQ)2) and v = lcm(2/i2, d{g)2) 
Ifg > 0, then S{g,t) = A{g,t){l + E2{vo)U,). 
Ifh is odd, then5{g,t) = A{g,t){l + E2{v)Ui). 
Ifg < 0, 2\h and 2\t, we have 6{g,t) = A{g,t). 

Next assume g < 0, 2\{h,t). 
If h2 = 2 and 8\d{go), then 



M((,,()(i + ni) if 4|«2. 



In the remaining cases we have 

(A{g,t)/2 if 2t2|/i2; 



5(9, t) 



Aig,t)/3 ift2 = /i2; 

A{g,t){l-lUi) ift2 = 2/i2; 



[A{g,t){l + Ui) if4/i2|t2. 

Corollary 1. (GRH). Let g > 1 be a square free integer. Then 

6{g,t) = {1 + E2{lcm{2,d{g)2))B{g,t))E{t). (7) 

Proof. We have A{g,t) = S{l,t,l) = E(t) (see the remark following Lemma 
H]). Furthermore, if 2\g and 2 \ t, then Hi = —B{g,t) and Hi = B{g,t) otherwise. 
Since E2{lcm{2,d{g)2)) = if g\2 and 2 f t, we infer that E2(lcm(2, d(^)2))ni = 
E2{lcm{2,d{g)2))B{g,t). Now invoke the theorem. □ 

Corollary 2. (GRH). Ift is odd, then 

%,t) = A(<?,t)(l-i(l-(-l)™)ni). 
Remark. On putting t = 1 one obtains the classical result of Hooley 

Proof of Theorem [T] On distinguishing cases according to the value of d{g)2, 
Corollary [T] yields Table 1. From Table 1 one easily reads off that if 6{g, t) = 0, then 
2 \ t, g = l(mod 4) and g\t. In this case we have {g/p) = 1 for the primes p \ g 
with p = l(mod t) by the law of quadratic reciprocity and hence Ng^t is finite and 
so6{g,t) = 0. ' D 

The proof of Theorem |3] will be given in Section HJ It will make use of properties of 
Wagstaff sums that will be established in the next section. 

3. Bringing the Wagstaff sums in Euler product form 
Recall the definition of the Wagstaff sum 

S[h,t,m) = > -. — — . 

^ ' ' ^ ^ nt^int) 

m I nt 
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A trivial observation is that if the divisibihty condition forces n to be non-square 
free, then /i(n) = and hence S{h,t,m) = 0. This happens for example if m2 f 2^2 
(cf. Lemma H]). 

In case m = 1 it is easily written as an Euler product (here we use that /i and (p 
are multiplicative functions). 

Lemma 1. 1). We have 



p\t, hp\tp p\t ' 

In 'particular, S{h, t,l) = iff 2\h and 2\t. 
2). If2\h and2\t, then 

p\t, hp\tp ^ p\2t ' 

Proof. 1) We have 

^''^ tifit)^ nip{nt){t,h) tip{t)\i-\ pip{pt){t,h). 

where we used that the sum S{h, t, 1) is absolutely convergent and the fact that the 
argument in the second sum is a multiplicative function in n. The contribution of 
the primes dividing t to this product is 

H n u-i) n i^-^^-^ n a-^). 

p\t, ptp\hp p\t, hp\tp p\t, hp\tp 

where we used that f{t)/t = Y[p\ti'^ ~ ^/p)- P\t: then 

{pt,h)ip{t) {p,h) 



and part 1 follows. 
2) We have 



pip{pt){t,h) p(p-l)' 

fi{n){nt,h) fi{n){nt, h) 



S(h,t,2) = Tt^^^^^^ = -\2 ^ ^ . 

^-^ ntifiint) ^-^ ntifiint) 

2\n ' 2\n ' 

The latter sum has the same Euler product as S{h,t,l), but with the factor for 
p = 2 omitted. □ 

Remark. The above lemma and the definition of the Artin constant shows that 
E{t) = S{l,t, 1) and A = S{1, 1, 1). 

Write M = m/ (m, t) and H = h/ {Mt, h). Then we have [13], Lemma 2.1] 

Sih,t,m) = f^iM)iMt,h)Eit) n -T^n 2 ^ , U "^T 11 f ^ ~ • 

q\(M,t) ^ l\M ^ ^ q\(t,H) ^ q\H ^ ^ 

q\t q\M q\Mt 
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The parameter H can be avoided as the formula can be rewritten as 

^^^^^f^n^ n n tf-i)n.(.-)- w 

q\mth 91'. qtq\hq q\t, hq\tq q\h 

mq\tq mq\tq q\mt 

(In order to see this it is helpful to consider the cases mq\tq, that is Mq = 1, 
and qtq\mq, that is q\M, separately.) These formulae relate S{h,t,m) to S{l,t, 1) 
(= E{t)), respectively to 1, 1) (= ^4), however, as we will show, expressions 
simplify considerably if we relate S{h,t,m) to S{h,t, 1). We start by showing how 
to remove odd prime factors from m. 



Lemma 2. Suppose that p f 2m. Then 
S{h, t, mp) — 



-5(/.,t,m)/(^-l) iip\t- 
S{h,t,m) if 



Proof. If p\t the summation condition mp\nt in the definition of S{h,t,mp) is 
equivalent with m\nt, that is we have S{h,t,mp) — S{h,t,m). 
Next assume that p\ t. We have 

S(h t mp) = V ^(^)(^^' ^) = fi{pn){pnt, h) ^ {p, h) ^i{n){nt, h) 

p\n p\n 

On noting that the latter sum can be written as S{h,t,m) — S{h,t,mp), the proof 
is then completed. □ 

Lemma 3. Suppose that we are not in the case where h is even and t is odd. We 

have 



S{h,t,2t2) 



-S{h,t, l)/3 if Icm(2,/i2)|t2; 
-S{h,t,l) if lcm(2, /i2) t ^2- 



Proof. We can write 

Sih t 2t ) - V ^ 1 ii{n){2nt,h) ^ ii{n){nt,h) 

where e is easily determined (and e ^ —1). Since the latter sum is equal to S{h, t, 1) — 
S{h, t, 2t2), we then infer that S{h, t, 2^2) = j^S{h, t, 1). Working out the remaining 
details is left to the reader. □ 

Lemma 4. Let m be an integer, having square free odd part. Let h and t be integers, 
with the requirement that t be even in case h is even. Then 

S{h,t,m)^ S{h,t,l)E^{m2) ~^ 



p(p-i) _ ^ ' 

P\m,p\2t (p,h) 
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where 

if m2\t2 

if 1712 = 2)^2 and lcm(2, /i2)|^2 
if 1712 = 2/f:2 and lcm(2, h2) \ ^2 
if m2 t 2^2, 
E2{m2), where E2{m2) is given by (jl]). 

Proof. By Lemma [T] the conditions imposed on h and t imply that S{h,t, 1) ^ 0. 
By Lemma 12] it suffices to show that S{h,t,m2) = S{h,t,l)Ei{m2). If ^21^2? then 
no divisibihty condition on n is imposed in the definition of S{h, t, m2) and so we 
obtain S{h,t,m2) = S{h,t,l) and hence Ei{m2) = 1. In case m2 = 2^2 we invoke 
Lemma |3l If m2 \ 2t2, then the summation condition m\nt implies 4|n and hence 
/i(n) = and so S{h,t,m2) = and hence Ei{m2) = 0. 

The final claim follows on noting that if 2h2\m2 and m2 = 2^2, then /i2|^2 and 
hence lcm(2, /i2) f iff 2 f ^2- □ 

4. Proof of Theorem E] 

The idea of the proof is to express S{g,t) in terms of S{h,t, 1), except in case 
g < 0, 2\h and 2 ft, when S{h, t, 1) = 0, in which case we express S{g, t) in terms of 
S{h,t,2). These two Wagstaff sums are then related to A{g,t) using the following 
lemma. Note that it shows that the dependence of A{g,t) on g is weak, as only h 
and the sign of g matter. 

Lemma 5. We have 

j^^g ^■^^i-S{h,t,2)/2 if (7<0, 2\h, 2\t- 
\S{h,t,l) otherwise. 

Proof. Note that li g < {] and 2|/i, then F2 = Q{i) and [F2 : Q] = 2. In the 
remaining cases we have [Fp : Q] = p{p — l)/{p,h). On invoking Lemma [T] the proof 
is then completed. □ 
Proof of Theorem^ Equation ([5]) follows by Lemma |5] and (|8]). We will use a few 
times, cf. the proof of Lemma [5l that 

n [ F : Ql - 1 = n £{P^_i- 

p\d{go), p>^t^ P p\d{go), p\2t ip,h) 

Assume GRH. 
The case g > 0. 

By Theorem [2] we have S{g,t) = S{h,t,l) + S{h,t,m), with m = lcm(2/i2, c?((7o)). 
First assume that 2\h and 2 \ t. Then, by Lemmas [1] and [5l we have S{h,t,l) = 
A{g,t) = and we need to show that S{g,t) = 0. Since S{h,t, 1) = it remains to 
show that S{h,t,m) = 0. Since for the n in the summation we have 4\2h2\n, this is 
clear. Next assume we are in the remaining case, that is either h is odd, or 2\{h,t). 
Then S{h,t,l) = A{g,t) by Lemma El Note that m2 = vq. By Lemma |4] we then 
find that 6{g,t) = S{h,t, 1)(1 + Ei{vo)Ui) = A{g,t){l + E2{vo)Ui), where we have 
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used that 2/i2|fo- 
The case h is odd. 

If (7 > then V = Vq and we are done, so assume that (? < 0. The formula for 
m in Theorem [2] can be rewritten as lcm(2, [(/(f/)!), and one finds that 6{g,t) = 
S{h,t,l) + S{h, t,lcm{2, \d{g)\)). This is the same formula as in case g > and 
2\ h, but with d{go) replaced by |c?((7)|. On noting that the odd part of d{go) equals 
the odd part of d{g), the result then follows. 
The case g<0, 2\t and 2\h. 

We have S{h,t,l) = S{h,t,m) = S{h,t,2h2) = and hence 6{g,t) = -S{h,t,2)/2 
by ([3D- Now invoke LemmaOto obtain 6{g,t) = A{g,t). 
The case g<0 and 2\{h,t). 

Note that 2|m and S{h,t,l) = A{g,t). By Lemma H] we infer that S{h,t,2) = 
S{h, t, 1) and S{h, t, 2/2.2) = S{h, t, 1)^2(2/12), where 

'1 if2/i2|t2; 
^2(2/12) = { -1/3 if h2 = t2] 
,0 if/l2t^2. 

Note that 

'1 if4|t2; 



^2(4) 



-1/3 ift2 = 2. 



If /i2 = 2 and 8|(i((7o); then by Theorem [2] we have m = 2gi, which can be rewritten 
as m = d{gQ)/2 (thus m2 = 4) and so 

t) = Sih, t, 1) - + ^^Mil + s{h, t, ^) 



5(/.,t,i)(^ + ^ + i?2(4)ni). 



where we used that, by Lemma HI S{h,t,d{go)/2) = S{h,t,l)E2{4:)Ili. Using that 
S{h,t, 1) = A{g,t) and the formula for -^2(4), we then arrive at ([6]). 
In the remaining case, m = lcm(4/i2, d{gQ)). Note that m2 = 4:h2 and 

1 if 4/i2|t2; 

£^2(4/12) = { -1/3 if 2h2 = t2] 
ii2h2\t2. 

We find that 

t) = ^(/., t, 1) - ^-^^ + ^Jhl^ + s{h, t, m) 

= 5(/.,M)(^ + ^^ + i?.(4/.2)n, 

Using that S{h, t, 1) = A{g, t) and the formulae for E2{2h2) and £^2(4/12) given above, 
the proof is then completed. □ 
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5. Vanishing of S{g,t) 

The aim of this section is to give a new proof of Theorem H] (due to Lenstra |3], 
who stated it without proof). The first pubhshed proof was given by Moree in [6]. 
He introduced a function Wg^t{p) G {0, 1, 2} for which he proved (see [6], for a rather 
easier reproof see [7]) under GRH that 

AT f \ (h +\ ( M(P-^)/^) , ^/a;logloga;\ 

Ng^tix) = {h,t) > Wg^tiP) + 0[—-^ ). 

, , . p - 1 V log X ' 

p<x, p=l(mod t) 

This function Wg^tip) has the property that, under GRH, Wg^tip) = for all primes 
p sufficiently large iff Ng^t is finite. Since the definition of Wg^t{p) involves nothing 
more than the Legendre symbol, it is then not difficult to arrive at the cases 1-6. 
E.g. in case 1 g is a square modulo p, and thus 2\t, contradicting 2 \ t. Likewise for 
the other 5 cases the obstructions can be written down (it turns out rg{p)2 7^ ^2 in 
each case). For the complete list of obstructions we refer to Moree [6l pp. 170-171]. 

Regarding the six vanishing cases Wagstaff 1^13, p. 143] wrote: 'It is easy to verify 
directly that our expression for S{g,t) vanishes in each of Lenstra's cases, but it is 
tedious to check that these are the only cases in which it vanishes'. We will show 
that once Wagstaff's result is brought into Euler product form, as done in Theorem 
[31 it is straightforward to establish Theorem HI A more conceptual, shorter and 
elegant (but less elementary) proof of Theorem |4] will appear in [5]. 

Theorem 4. (GRH). The set Ng^t is finite iff 6{g,t) = iff we are in one of the 
following six (mutually exclusive) cases: 

1) 2\t, d{g)\t. 

2) g>Q, 2h2\t2, ^\t, ^\h, d{-'ig^)\t. 

3) <7 < 0, /i2 = 1, t2 = 2, 3 1 1, ?.\h, d{?>g^)\t. 
A)g<^, h2 = 2, t2 = 2, d{2go)\2t. 

5) ^ < 0, /i2 = 2, t2 = 4, 3 1 1, 3\h, d{-Qgo)\t. 

6) g <0, 4/i2|t2, 3 ft, 3|/i, d{-3go)\t. 

Example. (GRH). If (7 > 1 is square free, then case 1 is the only one to take into 
account and we find S{g,t) = iff 2 f t, d{g)\t, that is iff 2 f t, g\t, g = l(mod 4). 



Table 2: Examples of pairs {g,t) satisfying cases 1-6 





1 


2 


3 


4 


5 


6 




(5,5) 


(3^4) 


(-l5^lo) 


(-6^6) 


(-6^4) 


(-3^4) 



Proof of Theorem |H If one of 1-6 is satisfied, then A^^^^ is finite. This can be 
shown by elementary arguments only involving quadratic reciprocity (see Moree [6], 
pp. 170-171]). It is thus enough to show that S{g,t) = iff one of the six cases is 
satisfied. For the proof we will split up case 6 into two subcases: 
6a) ^ < 0, 2|/i2, 4/i2|t2, 3 ft, 3\h, d{3go)\t. 
6b) g <0,h2 = l, 4|t2, 3 ft, 3\h, d{3go)\t. 

(For our proof it is more natural to require d{3go)\t, which, since 4\t, is equivalent 
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with d{—3go)\t.) Let us denote by d*{go) the odd part of the discriminant of go, that 
is d*{go) = d{go)/d{go)2- Note that 

fl iid*{go)\t- 
U,= l-1 if 3\d{go), d*{go)\3t, Sft, 3\h; (9) 

Ig(— 1,1) otherwise. 

The case 2 \t. 

If 2\h one has 5{g,t) = iS g > 0, that is iff d{g)\t. 

If 2 \ h, then A{g,t) ^ and we have b{g,t) = iff E2(lcm(2, (i(5()2)) = -1 and 

III = 1, that is iff \cm{2,d{g)2) = 2 and d*{g)\t, that is iff d{g)\t. 

Thus from now on we may assume that 2\t. This ensures that A{g,t) ^ 0. 

The case g > and 2\t. 

Now the possibihty E2{m2) = —1 cannot occur and thus 6{g,t) = iff £2(7712) = 1 
and Hi = —1. The latter two conditions are both satisfied iff lcm(2/;,2, '^(fi'o)2)|^2, 
3\d{go), d*{go)\3t, 3 ft, 3\h. These conditions can be reformulated as 2/i2|t2, 
3K^o), d{go)\3t, 3 t t and 3\h. Since 3 f t, 3\d{go), d{go)\3t iff d{-3go)\t, 3 \ t, 
we are done. 

Thus if (7 > or 2 1 1, then 5{g, t) = iff we are in case 1 or in case 2. It remains 
to consider the case where g < and 2\t. 
The case g <0, 2\t, 2\h. 

Here we have S{g,t) = iff E2{v) = 1 and Hi = —1. Note that E2{v) = 1 means 
that we require lcm(2, ^((7)2)1^2- 

If t2 = 2, then lcm(2, d{g)2)\t2 and Hi = —1 iff we are in case 3. 
If 4|t2, then lcm(2, d{g)2)\t2 and Hi = —1 iff we are in case 6b. 
The case g <0, 2\{h,t). 

We have 6{g, t) = iff we are in one of the following three cases: 

A) h2 = 2,t2 = 2, 8\d{go), Hi = 1; 

B) h2 = 2, t2=4, 8|rf(^o), ni = -l; 

C) 2|/l2, 4/l2|t2, Hi = -1. 

It is easily checked that these are merely cases 4, 5 and 6a in different guises. 

To sum up, we have shown that 6{g, t) = iff we are in one of the cases 1,2,3,4,5,6a 
or 6b. Note that the six cases are mutually exclusive. □ 

We now propose a conjecture on 6{g,t) for arbitrary rational g. It generalizes 
Conjecture 1. 

Conjecture 2. The set Ng t has a natural density 6{g, t) that is given as in Theorem 
and is a rational multiple of the Artin constant A. The set Ng^t is finite iffS{g, t) = 
iff we are in one of the six cases of Theorem^ 

On combining Theorem [3] and Theorem |4] we deduce that Conjecture [2] holds true 
on GRH. 

Theorem 5. Conjecture\^is true under GRH. 
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6. Near-primitive roots density through character sum averages 

Lenstra, Moree and Stevenhagen ^ show that for a large class of Artin-type 
problems the set of primes has a natural density S that is given by 

6=il + llEp)llAp, (10) 
p p 

where Yl^ Ap is the 'generic answer' to the density problem (e.g. A in the original 
Artin problem) and 1 + Yip Ep a correction factor. For finitely many primes p one 
has Ep 1 and further — 1 < ii^p < 1 as £'p is a (real) character sum average over 
a finite set (and hence the correction factor is a rational number). In particular, it 
is rather easy in this set-up to determine when 6 = 0. The character sum method 
makes use of the theory of radical entanglement as developped by Lenstra |1] 

For the near-primitive root problem the method leads rather immediately to the 
formula S{g,t) = A{g,t){l + -^2-^1) in case 5^ > 0. The only harder part is the 
determination of E2. For the details the reader is referred to [5]. 

Indeed, the great advance of the newer method is that it very directly leads to a 
formula for the density in Euler product form. The classical method leads to infinite 
sums involving the Mobius function and nearly multiplicative functions (in our case 
Wagstaff's result (Theorem [2]). It then requires rather cumbersome manipulations 
to arrive at a density in Euler product form. Indeed, inspired by the predicted 
result (fTOl) the author attempted (and managed) to bring Wagstaff's result in Euler 
product form. 

The analogue of Theorem [3] obtained in this approach. Theorem 6.4 of [^Sj, looks 
slightly different from Theorem El However, on noting that S2 as defined in Theorem 
6.4 is merely the 2-part of m as defined in Wagstaff's result Theorem [2], it is not 
difficult to show that both methods give rise to the same Euler products for the 
density. By allowing to be negative in case h is odd and g < 0, the above 6 cases 
where vanishing occurs can be reduced to 5 cases (see Corollary 6.5 of [5]). 



7. An application 

Let $„(x) denote the n-th cyclotomic polynomial. Let 5* be the set of primes p such 
that if f{x) is any irreducible factor of ^p{x) over F2, then f[x) does not divide any 
trinomial. Over F2, ^p{x) factors into r2{p) irreducible polynomials. Let 

Si = {{p > 2 : 2\ r2{p)}} U {p > 2 : 2 < r^ip) < 16})\{3, 7, 31, 73}. 

Theorem 6. We have Si C S . The set Si contains the primes p > 3 such that 
p = ±3(mod 8). On GRH the set S\ has density 

1 ^1323100229 ^ 

^ ' 2 1099324800 ^ ' 

Proof. The set > 2 : 2 f r2(p)}} equals the set of primes p such that (^) = — 1, 
that is the set of primes p such that p = ±3(mod 8). This set has density 1/2. We 
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thus find, on consulting Table 1, that 
S{Si) = 1+ ^(2, 



2 ' / ^ --\-iJ) 

2<j<16 

= \ + + ^ + ll + ^) + ^(6)(1 + ^) + ^(10) + i5(14), 

which yields f|TT]) on invoking the definition ([1]) of -E(t). That Si C S* is a conse- 
quence of the work of Golomb and Lee [1] . □ 

Acknowledgment. Given expressions like ([H]), niy intuition was that expressing 5{g, t) 
in Euler product form would lead to very unpleasant formulae and thus I never at- 
tempted this. Discussions with Peter Stevenhagen, considering the near-primitive 
root problem by a much more algebraic method, strongly suggested easier expres- 
sions for 5{g, t) than expected. This led me to try to bring Wagstaff's result in Euler 
product form, also with the aim of verifying the results found by the character sum 
method (alluded to in Section [6]). 

The memorable discussions with Stevenhagen are gratefully acknowledged. Fi- 
nally, I like to thank Carl Pomerance for helpful e-mail correspondence. 
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